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Abstract 

Several  3-D,  viscous,  turbulent  Navier- 
Stokes  computations  have  been  performed 
for  a  missile  equipped  with  and  without 
divert  jet  thruster  for  three  different  wing 
planforms  (having  a  fixed  tail  configuration) 
at  a  nominal  flow  Mach  number  of  3.94, 
angles-of-attack  ranging  from  2  to  25  deg 
and  jet  thrust  ratios  of  one  and  four.  Results 
are  presented  to  show  that  the  normal  force 
and  pitching  moment  coefficients  for  aU 
computed  cases  are  predicted  within  5 
percent  of  the  wind-tunnel  data.  Synthesis  of 
aU  the  results  show  low  amplification  factor 
for  all  windward  jet  thruster  cases.  For  this 
case  the  upstream  favorable  pressure  zone 
created  by  the  windward  jet  is  insufficient  to 
compensate  for  the  massive  unfavorable 
pressure  loss  on  the  windward  wings/tails 
due  to  the  jet  blockage  and  jet  wrap-around 
effects. 


Nomenclature 

AF  Amplification  Factor, 

1  +  (CNje,  -  CN„o.je,)/(T/q«S) 

Alpha  Angle-of-Attack  (deg) 

*Senior  Development  Engineer 
System  Design  Laboratory, 
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CN  Normal  Force  Coefficient, 

CN=.(N/q.S) 

dp  Pressure  Differential, 

dp  =  (Pjet  -  Pno-jet)/Gamma  •  Rnf) 
Gamma  Ratio  of  Specific  Heats 
M  Freestream  Mach  Number 

N  Normal  Force  (N,  lb) 

P  Pressure  (N/m^,  Ib/ft^) 

Phi  Azimuth  Angle  (deg) 

q  Dynamic  Pressure,  q  =  1/2  pv^ 

S  Missile  Cross-Sectional  Area 

T  Jet  Thrust  (N,  lb) 

V  Velocity  (m/sec,  ft/sec) 

Greek 


cc  Angle-of-Attack  (deg) 

p  Density 

(j)  Azimuth  Angle  (deg) 

Subscript 

inf  Freestream  condition 

jet  Condition  with  jet 

no-jet  Condition  with  no  jet 

Introduction 

Rapid  airframe  response  time  is  critical 
during  the  homing  phase  of  interceptor 
missiles.  Surface  mounted,  fast-reacting  jet 
thrusters  offer  an  attractive  alternative  to 
conventional  aerodynMiic  surface  control  for 
improved  missile  agility  and  maneuverability. 
Performance  enhancements  at  low  speed  and 
high  altitude  (where  dynamic  pressure  is 
low)  are  additional  advantages  of  the  surface 
reaction  jets.  It  has  been  found  that  under 
some  missile  orientations  and  flow  condi¬ 
tions,  the  surface  mounted  lateral  jet  leads  to 
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thrust  amplification  due  to  a  high  surface 
pressure  region  that  forms  ahead  of  the  jet 
on  the  missile  surface.  Effective  exploitation 
of  the  thrust  amplification  can  lead  to 
improved  missile  performance.  However,  it 
has  also  been  observed  that  at  certain  other 
flow  conditions  and  orientations,  the  reaction 
jets  can  produce  negative  effects  (i.e.,  thrust 
deamplification).  An  understanding  of  the 
controlling  factors  that  produce  high  amplifi¬ 
cation  as  well  as  those  that  produce  low 
amplification  are  critical  in  developing  a 
credible  design  basis  for  optimal  missile 
performance. 

The  problem  of  the  lateral  jet  interaction 
with  the  external  flow,  under  conditions  of 
varying  flight  Mach  numbers,  angles-of- 
attack  and  jet  orientations,  is  extremely 
complex  in  nature  and  has  been  studied 
experimentally  and  analytically  for  many 
years^’^.  More  recently.  Computational  Fluid 
Dynamics  (CFD)  studies^'’  have  been 
performed  to  understand  these  effects.  Our 
approach  is  to  judiciously  combine  wind- 
tunnel  testing  and  CFD  simulation  in  an 
effort  to  evolve  a  validated  design  and 
analysis  tool  that  can  S5mthesize  the  physical 
complexity  of  the  flow  and  identify  its  key 
controlling  parameters.  This  paper  deals 
with  the  overall  CFD  validations  with  and 
without  lateral  jet  thrusters  and  subsequently 
outlines  the  utility  of  the  approach  for  later^ 
jet  effectiveness  studies  in  relation  to  a 
missile  design. 

This  paper  is  subdivided  into  several 
sections.  The  next  section  briefly  outlines 
the  previous  work  in  this  area  using  CFD 
approaches.  The  details  of  the  computa¬ 
tional  methodology,  geometry,  grid  related 
issues  and  boundary  conditions  for  the  CFD 
application  are  discussed  later.  CFD  valida¬ 
tion  studies  for  a  generic  class  of  missile 
airframes  with  and  without  lateral  thrusters 


are  then  compared  with  the  wind-tunnel  data 
for  the  normal  forces  and  moments.  A  later 
section  presents  the  overall  summary  and 
conclusions  for  the  paper. 

Background 

The  topic  of  jet  interaction  with  an 
external  supersonic  flow  dates  back  to  the 
mid- 1960^’^,  when  a  large  number  of  generic 
experimental  data  were  generated  and  related 
correlation  techniques  were  developed. 
Emergence  of  hypersonic  interceptors, 
maturity  of  CFD  and  the  advent  of  super¬ 
computers  revived  these  activities  in  the  late 
1980’s^'^.  Several  investigators  have  per¬ 
formed  CHD  studies  for  the  fundamental 
problem  of  jet  interaction  in  relation  to 
adaptive  gridding®,  turbulence  models^,  grid 
refinements^®,  and  the  impact  of  artificial 
viscosity*.  These  studies  range  from  Euler^^ 
to  Navier-Stokes  computations^^.  Of  these 
studies,  particular  reference  is  made  to  the 
studies  reported  by  Dash  etal.^^  and  York 
et.al.^  because  much  of  the  current  CFD 
effort  is  derived  from  their  mature  technical 
expertise  in  this  area.  Further  details  of  the 
methodology  and  related  research  work  can 
be  obtained  from  the  references  cited. 

While  a  vast  number  of  numerical  studies 
have  been  performed  using  controlled  jet 
interaction  studies  for  methodology  devel¬ 
opment,  efforts  to  simulate  missile  surfaces 
have  been  rather  limited.  More  recently, 
Chan  et.al.^*  performed  a  series  of  studies 
that  lead  to  the  simulation  of  a  full  missile 
surface  with  control  surfaces  and  jet  interac¬ 
tion.  Qin  and  Foster^"*  also  performed  similar 
studies  using  a  Navier-Stokes  approach  for 
an  inclined  jet  on  an  ogive/cylinder  body. 
These  results  depict  the  remarkable  flow 
details  obtained  using  CFD  approaches. 
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which  ultimately  result  in  making  judicious 
choices  for  flight  vehicles  design  and  further 
wind-tunnel  testing.  Srivastava^^  performed 
full  Navier-Stokes  studies  for  generic  missile 
bodies  with/without  leeward  and  windward 
jets  but  without  wing  or  tail  panels.  Com¬ 
parisons  with  the  wind-tunnel  data,  however, 
was  not  direct  because  the  tests  were  con¬ 
ducted  with  tail  panels  while  computations 
were  performed  without  tail.  Removal  of  the 
tail  load  from  the  wind-tunnel  data  by  an 
approximate  method  introduced  uncertainties 
that  were  not  fully  quantified.  This  defi¬ 
ciency  in  the  CFD  model  has  been  eliminated 
in  the  current  paper,  which  shows  direct 
comparison  of  the  CFD  predictions  with  the 
wind-tunnel  data  by  modeling  all  geometrical 
aspects  of  the  wind-tunnel  missile  geometry. 

Computational  Methodology 

PARCH^,  which  is  a  fuU  Navier-Stokes 
(FNS)  code  with  plume/missile  airframe 
steady-flow  predictive  capability,  is  being 
used  for  our  current  studies.  PARCH  code 
utilizes  formulations  based  on  the  NASA/ 
Ames  ARC  aerodynamic  code  and  the 
AEDC  propulsive  extension,  PARC.  This 
code  is  particularly  suited  for  missile  surfaces 
due  to  its  grid  patching  capability  which  is 
useful  for  treating  embedded  surfaces  in  a 
flow  field.  Patching,  that  is  accomplished  in 
mapped  computational  coordinates,  is 
automatically  constructed  from  boundary 
inputs.  Boundary  conditions  are  applied 
along  the  outer  computational  boundaries 
and  relevant  embedded  surfaces.  The  code 
utilizes  diagonalized  Beam-Warming  numer¬ 
ics  with  matrix- split  finite  rate  chemistry. 
Several  versions  of  the  K-E  turbulence 
model  are  available  in  the  code  that  were 
specifically  developed  for  jet  interaction  and 


propulsive  studies.  We  are  currently  using 
the  capped  low  Reynolds  number  formula¬ 
tion  of  Qhien’s  K-E  model^^  for  the  current 
simulations.  Further  details  of  the  code 
capability  can  be  found  in  Reference  6. 

Typical  boundary  approaches  for  the 
current  application  (supersonic  flows)  are 
specified  supersonic  freestream  conditions  at 
the  inlet  and  outer  boundaries.  Extrapolation 
procedures  are  employed  at  the  exit  bound¬ 
ary.  Surface  conditions  are  appropriate  to 
viscous  flows  with  adiabatic  wall  condition. 

Surface  jet  boundary  condition  is  the 
specified  jet  nozzle  supersonic  exit  condi¬ 
tions.  The  circular  area  of  the  jet  in  the 
wind-tuimel  test  is  approximated  by  a  square 
aperture  in  the  CFD  simulation. 

Fig.  1  shows  a  sketch  of  the  wind-tunnel 
models  recently  tested  by  Raytheon  in  McAir 
Polysonic  Wind-Tunnel  with  three  different 
wing  planforms.  The  wing  in  Fig.  la  is  the 
baseline  configuration  which  is  15.734  in. 
long  and  its  leading  edge  starts  at  12.134  in. 
from  the  missile  nose.  The  jet  thruster  is 
located  at  14  in.  from  the  nose.  The  jet  area 
is  0.1764  in.^  circular  hole.  The  second 
configuration,  shown  in  Fig.  lb,  is  a  forward 
wing  configuration  whose  length  remains  the 
same  but  its  leading  edge  is  moved  forward 
to  10.301  in.,  from  the  missile  nose.  The 
third  configuration  is  a  short  wing,  length 
7.67  in.  with  leading  edge  starting  at  10.301 
in..  For  aU  these  configurations,  tail  panel 
geometry  as  well  as  the  axial  jet  location  are 
the  same.  Missile  is  rolled  in  the  tunnel  to 
achieve  windward  and  leeward  jet  orienta¬ 
tions.  The  fuU  CFD  simulation  geometry  is 
shown  in  Figs.  2-4  for  the  three  configura¬ 
tions.  CFD  simulations  utilize  the  rotation  of 
inlet  velocity  vector  to  simulate  leeward  and 
windward  jet  orientations.  These  simulations 
were  performed  on  three  different  grids 
corresponding  to  each  configuration.  The 
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axial  distribution  of  the  grid  was  the  only 
component  that  was  modified  keeping  the 
total  number  of  grids  the  same,  i.e., 
231*51*69,  axial,  normal  and  circumferential 
in  that  order.  Grid  patching  procedure  was 
used  to  apply  the  relevant  boundary  condi¬ 
tions  on  the  surfaces.  AU  angles-of-attack 
for  a  given  configuration  were  simulated  on 
a  single  grid  consistent  with  the  minimum 
desired  Mach  number  and  maximum  desired 
angle-of-attack.  This  allows  us  to  minimize 
on  the  grid  effort.  The  grids  were  generated 
through  ‘GRIDGEN”,  with  geometry  models 
developed  within  ‘GRIDGEN”^^ 

•  The  jet  area  was  resolved  by  using  a  10  x 
10  grid  on  the  missile  surface.  The  boundary 
conditions  on  this  jet  area  grid  were  specified 
by  inputting  the  exit  conditions  of  a  nozzle 
analysis  using  an  iterative  procedure  consist¬ 
ing  of  the  nozzle  exit  area,  isentropic  condi¬ 
tions  and  the  specified  jet  thrust.  Local  time 
stepping  procedures  with  conventional  dissi¬ 
pation  factors  were  employed  to  achieve 
convergence  defined  by  at  least  three  orders 
of  magnitude  reduction  hi  the  maximum 
residual.  Converged  solutions  were  checked 
with  three  and  four  orders  of  magnitude 
reduction  in  the  maximum  residual  for 
selected  cases.  The  minimum  grid  spacing 
near  the  missile  surface  was  assigned  by 
choosing  a  y"^  value  of  5  for  all  computations. 
Variation  of  this  value  did  not  produce  any 
significant  changes  in  the  predicted  results. 
Post-analysis  of  the  computed  cases  later 
suggested  that  the  external  flow  is  primarily 
supersonic  and  the  vortical  flow  nature 
around  the  missile  geometry  is  primarily 
governed  by  the  inviscid  phenomena. 
Selected  test  cases  with  inviscid  CFD 
simulations  showed  that  force  and  moment 
coefficients  can  be  well  predicted  using 
inviscid  models. 


Comparison  With  Wind-Tunnel  Tests 

Our  current  effort  was  confined  to  a 
nominal  flow  Mach  number  of  3.94  with 
angles-of-attack  varying  fi'om  2  to  25  deg. 
Exact  wind-tunnel  conditions  were  used  for 
each  CFD  simulation.  CFD  efforts  were 
restricted  to  leeward  and  windward  jet 
orientations,  which  encompass  the  most 
severe  flow  cases.  This  also  allows  us  to 
exploit  symmetry  in  CFD  simulation.  Jet 
thrusts  were  simulated  at  the  wind-tunnel 
test  conditions  of  175  and  50  lb  (nominal 
cases). 

Baseline  Wings 

Windward  Jets 

Fig.  5a  shows  the  computed  surface 
pressure  contours  for  two  cases.  The  top 
figure  is  with  a  windward  jet,  jet  thrust  175 
lb,  at  a  flow  Mach  number  of  3.94  and  at  an 
angle-of-attack  of  20  deg.  For  this  case, 
notice  the  extent  of  high  pressure  forward  of 
the  jet  and  very  low  pressures  aft  of  the  jet 
on  the  windward  side.  It  is  quite  clear  fi:om 
this  figure  that  the  presence  of  the  jet-on  the 
surface  has  totally  wiped  out  the  aft  normal 
forces  created  by  the  windward  wings.  This 
is  more  clear  if  we  compare  the  lower  figure 
which  shows  jet-off  pressure  contours,  with 
the  jet-on  upper  figure.  The  table  on  this 
figure  shows  the  comparison  of  the  normal 
force  and  pitching  moment  coefficients  with 
the  wind-tunnel  data.  The  comparisons  are 
excellent  Overall  we  can  conclude  from  this 
figure  that: 

(1)  Windward  located  jet  thruster  cre¬ 
ates  a  favorable  high  pressure  zone 
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ahead  of  the  jet,  but  only  at  the  ex¬ 
pense  of  a  massive  loss  of  normal 
force  from  the  windward  wings. 

(2)  For  this  case  the  thrust  amplifica¬ 
tion  factor  (defined  as  1  +  aero- 
interactive  force/thrust,  aero- 
interactive  force  =  jet-on  normal 
force  on  the  airframe  -  jet-off  nor¬ 
mal  force  on  the  airframe)  is  0.56. 
Note,  however,  that  the  pitching 
moment  is  significantly  increased 
with  the  jet  thrusters.  Even  though 
the  latter  can  be  used  effectively  to 
our  advantage  for  vehicle  maneu¬ 
ver,  our  goal  is  to  also  enhance  the 
amplification  factor. 

(3)  Notice  also  from  these  figures  that 
a  jet  wrap-around  effect  is  seen  to 
create  an  adverse  effect  causing 
undesirable  normal  forces. 

(4)  The  effect  of  the  jet  gas  on  the  tail 
panels  is  also  massive,  resulting  in 
the  loss  of  tail  control  power  for 
this  jet  orientation 

These  effects  were  observed  in  our  previous 
computations'^  but  now  the  results  are  far 
more  credible  due  to  the  excellent  compari¬ 
sons  with  the  wind-tunnel  data. 

The  vortical  flow  pattern  of  the  missile 
with  and  without  the  windward  jet  can  be 
observed  from  Figs.  5b- 5f.  In  Fig.  5b,  the 
total  pressure  contours  for  several  axial 
locations  of  the  missile  without  the  divert 
thruster  is  shown.  Notice  the  early  forma¬ 
tion  of  the  body  vortices  on  the  leeward  side, 
forward  of  the  wings,  which  interact  with  the 
leeward  wings,  are  subsequently  modified 
and  later  interact  with  the  tail  panels. 
Additional  vortices  are  generated  on  the 
leeward  side  of  the  lower  wings  that  later 
follow  an  upward  path  as  seen  from  this 
figure.  Fig.  5c  shows  the  Mach  number 


contours  for  the  same  axial  locations  as  in 
Fig.  5b  showing  fully  supersonic  flows  for 
the  entire  flow  field.  For  these  flow  condi¬ 
tions,  the  inviscid  effects  dictate  the  vortical 
flow  structure. 

Figs.  5d-5f  show  similar  plots  with  a 
windward  divert  thruster.  Notice  from  these 
figure  that  the  jet  effects  result  in  massive 
loss  of  the  total  pressure,  first  on  the  wind¬ 
ward  side  but  gradually  wrapping  around  the 
missile  downstream  to  engulf  the  leeward 
side  of  the  missile  surfaces. 

Fig.  6  shows  similar  computations  for  a 
lower  jet  thrust  of  50  lb  at  identical  flow 
conditions.  Once  again,  first  observe  the 
excellent  comparison  with  the  wind-tunnel 
data.  For  this  case  the  thrust  amplification 
factor  is  -0.17,  indicating  that  the  thruster 
loses  its  effectiveness  except  for  the  gains  in 
the  moment  coefficient  Figs.  7a  and  7b 
show  the  previous  results  in  a  different 
format.  In  an  effort  to  enhance  the  interac¬ 
tive  effects,  a  pressure  differential  plot 
between  the  jet-on  and  jet-off  cases  are 
shown  for  high  and  low  thrust  cases.  Note 
from  these  figures  that  the  high  pressure 
zone  ahead  of  the  jet  wrap-around  effects  of 
the  jet  and  massive  loss  of  pressure  forces  on 
windward  wing  and  tail  panels  are  the 
controlling  physical  mechanisms. 

At  a  jet  thrust  of  175  lb,  the  effect  of 
angles-of-attack  are  shown  in  Figs.  8  and  9. 
Note  that  at  an  angle-of-attack  of  10.8  deg, 
comparisons  with  wind-tunnel  remain 
excellent  and  the  forward  high  pressure 
interactive  zone  seems  to  shrink  in  extent 
especially  in  the  front  end.  This  is  more 
visible  in  Fig.  9  for  angle-of-attack  of  4.02 
deg  with  excellent  comparisons  with  the 
wind-tunnel  data. 

Further  computations  were  made  for 
angles-of-attack  of  9.95  and  3.13  deg 
without  the  jet  thrusters  to  vahdate  the  CFD 
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results  at  low  angles-of- attack.  These  results 
are  shown  in  Figs.  10  and  11.  Predicted 
results  are  in  excellent  comparison  with  the 
experimental  data. 

Leeward  Jets 

Figs.  12-15  show  the  CFD  computations 
for  leeward  jet  orientation.  Two  thrust 
levels  were  computed,  175  and  50  lb.  The 
angle-of- attack  range  was  from  2.2  to  20  deg 
at  a  nominal  flow  Mach  number  of  3.94.  All 
the  computed  results,  as  shown  in  these 
figures,  show  excellent  comparisons  with  the 
wind-tunnel  data.  Fig.  12a  shows  two  views 
of  the  computed  results  for  a  jet  thrust  of 
175  lb  at  an  angle-of-attack  of  19  deg.  View 
I  shows  the  leeward  side  while  view  n  shows 
the  windward  side.  A  comparison  of  this 
figure  with  the  results  shown  in  Fig.  5  show 
the  following: 

(1)  For  a  leeward  jet  thruster,  there  is 
only  a  marginal  positive  pressure 
gain  ahead  of  the  jet  because  the  jet 
exhausts  in  a  vortical  flow  zone. 

(2)  The  windward  wings  and  missile 
body  pressure  levels  are  also  af¬ 
fected,  as  is  seen  if  we  compare 
Fig.  12a  view  n,  with  the  jet-off 
windward  side  in  Fig.  5. 

(3)  Comparison  with  the  experimental 
data  is  excellent 

Figs.  12b-12d  show  the  vortical  flow 
pattern  for  the  leeward  divert  thruster. 
Notice  from  the  total  pressure  contours  of 
Fig.  12b  that  the  effect  of  the  leeward  jet  is 
confined  between  the  two  leeward  wings. 
Comparison  of  this  figure  to  Fig.  5b  (without 
the  divert  thruster),  shows  that  the  windward 
flow  pattern  for  the  two  cases  are  very 
similar.  This  indicates  that  for  leeward  jets. 


jet  effects  are  confined  only  to  the  leeward 
side  of  the  missile  surfaces.  Fig.  12d  shows 
the  Mach  number  contours  for  this  case 
showing  predominantly  supersonic  flow 
everywhere  in  the  flow  field. 

Fig.  13  shows  a  similar  comparison  for  a 
leeward  jet  with  a  jet  thrust  of  50  lb  and  an 
angle-of-attack  of  20  deg.  The  table  shows 
the  excellent  comparison  of  the  computed 
results  with  the  experimental  data.  Notice 
also  from  this  figure  that  the  windward 
control  surface^  are  nearly  unaffected  by  the 
presence  of  the  leeward  jets.  Figs.  14  and  15 
show  that  at  lower  angles-of-attack,  9  and 
2.2  deg,  respectively,  the  positive  pressures 
ahead  of  the  leeward  jets  are  somewhat 
recovered  due  to  the  interaction  with  ap¬ 
proaching  streamlines  over  the  missile  body. 

Forward  Wings 

In  an  attempt  to  understand  the  effects  of 
the  jet  thruster  location  in  relation  to  the 
missile  wings,  tests  were  conducted  with 
wings  moved  forward  of  the  jet  location  by 
1.83  in.  The  objective  was  to  attempt  to 
enhance  the  positive  pressure  effects  forward 
of  the  jets,  for  windward  jet  orientation,  in 
an  effort  to  Enhance  amplification  factors. 
Fig.  16  shows  the  computations  with  a 
forward  wing,  with  an  without  windward 
jets,  at  an  angle-of-attack  of  25  deg  and  a 
nominal  flow  Mach  number  of  3.94.  Once 
again,  the  excellent  comparisons  between  the 
computed  results  and  wind-tunnel  data  is 
obvious  from  the  table  on  this  figure.  It  is 
also  observed  from  this  figure  that  with  a 
forward  wing,  the  positive  pressure  zone  is 
somewhat  larger  as  compared  to  Fig.  5  for 
the  baseline  wing  (note  that  the  angle-of- 
attack  is  lower  in  Fig.  5).  The  amplification 
factor  for  this  case  is  0.47.  This  happens 
because  the  loss  of  positive  pressure  on  the 
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windward  wings  due  to  the  jets  effects  still 
dominates  the  gain  ahead  of  the  jet.  This 
phenomena  persists  at  a  lower  jet  thrust  level 
of  102  lb  for  this  case  as  shown  in  Fig.  17. 

Short  Wings 

A  short  wing  (length  7.672  in.  as  com¬ 
pared  to  baseline  case  of  15.734)  missile  was 
also  tested  in  the  wind-tunnel  to  investigate 
jet  effects.  Fig.  18  shows  the  computed 
results  for  this  case  at  an  angle-of-attack  of 
20  deg  and  a  flow  Mach  number  of  3.94, 
with  and  without  jet  effects.  The  overall 
CFD  comparisons  for  normal  and  pitching 
moment  coefficients  are  stiU  excellent,  as 
shown  in  the  table  on  this  figure.  The 
amplification  factor  for  this  case  is  1.12 
which  compares  to  0.56  for  the  baseUne 
Aving  at  the  same  condition.  It  would  appear, 
at  least  superficially,  that  significant  gains  in 
amplification  factor  may  be  achieved  by 
careful  selection  of  wing  planform  and  the 
location  of  the  lateral  jet  relative  to  the  wing. 
The  same  behavior  is  observed  in  Fig.  19  at  a 
reduced  jet  thrust  level  of  50  lb.*  In  Fig.  20, 
the  pressure  differential  between  the  jet-on 
and  the  jet-off  conditions  are  shown  for  the 
short  wing.  This  figure  highlights  the  jet 
wrap-around  effect  and  loss  of  tail  load  with 
the  jet-on  conditions. 


Rolled  Baseline  Wings 

When  the  missile  is  roUed  such  that  one 
of  the  wings  is  vertically  downward,  the 
vortical  flow  pattern  is  likely  to  be  influ¬ 
enced.  The  computational  results  and 
comparison  with  experiments  for  this  missile 
orientation  is  discussed  next 

Fig.  21a  shows  the  computed  pressure 
contours  for  this  orientation  of  the  missile  at 


an  angle-of-attack  of  20  deg  and  a  flow 
Mach  number  of  3.94.  This  figure  also 
shows  the  comparison  of  the  experimental 
data  with  the  computed  force  and  moment 
coefficients.  Notice  that  the  experiments 
were  at  a  slightly  different  condition  that 
introduces  asymmetry  where  as  CFD  compu¬ 
tations  were  made  for  a  symmetric  case.  A 
somewhat  clear  view  of  the  pressure  con¬ 
tours  are  shown  in  Fig.  21b.  Notice  that  the 
horizontal  wings  are  highly  loaded  for  this 
orientation  of  the  missile.  Fig.  21c  shows 
the  total  pressure  contours  indicating  vortical 
flow  pattern  for  this  missile  orientation  while 
Fig.  2 Id  shows  the  Mach  number  contours. 
Notice  again  the  supersonic  nature  of  the 
flow  pattern  for  the  whole  flowfield. 

Overall  Results 

All  the  predicted  results  from  CFD  and 
wind-tunnel  tests  for  the  normal  force 
coefficients  and  the  pitching  moment  coeffi¬ 
cients  are  shown  in  Figs.  22  and  23,  respec¬ 
tively,  indicating  the  degree  of  deviation  of 
the  predicted  results  from  the  wind-tunnel 
data.  Notice  that  both  curves  show  very 
little  deviation  from  the  desired  45  deg  line. 
The  corresponding  table  of  data  is  shown  in 
Table.  I. 

Conclusions 

Eighteen  widely  varying  wind-tunnel  test 
cases,  encompassing  geometry  variations 
(baseline,  forward  and  short  wings),  angles- 
of-attack  (2  to  25  deg),  lateral  jet-on  with 
windward  and  leeward  orientations  and  jet- 
off,  were  simulated  using  CFD  to  show  that 
the  normal  force  and  moment  coefficients  for 
aU  cases  are  in  excellent  agreement  with  the 
data.  Fig.  22  for  the  normal  force  coefficient 
and  Fig.  23  for  the  pitching  moment  coeffi- 
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dent  sum  up  the  overall  accuracy  of  the 
predicted  results  for  all  the  cases  considered. 
These  validations  demonstrate  that  our 
current  CFD  methodology  can  provide 
accurate  predictions  of  the  force  and  moment 
coefficients  at  widely  varying  flow  condi¬ 
tions.  Additionally,  synthesis  of  the  CFD 
results  for  various  cases  highlight  the  physi¬ 
cal  mechanisms  that  inhibit  the  jet  amplifica¬ 
tion  factor  for  windward  jet  orientation. 
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TABLE  I  -  Table  Showing  Comparisons  of  CFD  and  Wind  Tunnel  Force 

and  Moment  Coefficients 


CASE# 

Mach 

Angle  of 

Jet  Thrust 

CN 

CN 

CM 

CM 

Number 

Attack-deg 

T  ( Lb.) 

CFD 

EXPT. 

CFD 

EXPT. 

1.  BW 

3.94 

20 

175 

8.96 

8.71 

22.56 

22.5 

2,  BW 

3.94 

20 

0 

7.3 

6.8 

5.96 

5.97 

3,  BW 

3.94 

20 

50 

6.91 

6.63 

11.96 

12.28 

4,  BW 

3.94 

10.8 

175 

6.54 

6.30 

16.03 

15.76 

5  ,  BW 

3.94 

4.02 

175 

4.86 

4.70 

13.48 

13.26 

6,  BW 

3.94 

9,95 

0 

2.83 

2.84 

2.71 

2.89 

7,  BW 

3.94 

3.13 

0 

0.73 

0.79 

0.91 

1.081 

8,  BW 

3.94 

19.0 

-175 

2.49 

2.35 

-6.09 

-  6.20 

9,  BW 

3.94 

20.0 

-50 

5.55 

5.46 

2.15 

2.24 

10,  BW 

3.94 

9 

-175 

-1.47 

-1.4 

-9.4 

-9.09 

11,  BW 

3.94 

2,2 

-175 

-3.39 

-3.3 

-11.59 

-11.37 

12,  FW 

3.94 

25.0 

175 

11.13 

10.77 

29.24 

29.70 

13,  FW 

3.94 

25.0 

0 

9.28 

9.17 

11.03 

11.66 

14,  FW 

3.94 

25.0 

102 

9.66 

9.40 

23.03 

23.52 

15,  SW 

3,94 

20.0 

175 

9.31 

8.97 

25.11 

25.32 

16,  SW 

3.94 

20.0 

50 

6.45 

6.02 

13.78 

13.68 

17,  SW 

3.94 

20.0 

0 

5.53 

5.15 

7.65 

7.49 

NOTE: 


BW  - 

Baseline  Wing 

FW  - 

Forward  Wing 

SW  - 

Short  Wing 

CN  - 

FN/(q*S) 

CM  - 

Fn  *(Xmc  -Xcp)/(q*S*X^,) 

Fn  - 

Normal  Force  ( Lbs) 

q 

Dynamic  head 

s  - 

Cylinder  Cross-  Sectional  Area 

Xrei  .  Reference  Length  (  Cylinder  Diameter) 

Xmc  .  Moment  Reference  distance 
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Fig.  2  -  Wind  Tunnel  Geometry  for  Computational  Fluid  Dynamics  Simulation 
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Fig.  3  -  Wind  Tunnel  Geometry  for  Computational  Fluid  Dynamics  Simulation 


UNCLASSIFIED 


sz 

o 

c. 


o 

CO 

o  ^ 
c5  .E  ■ 

.2>  ^  -D 
K  LU 

rl  '  cj) 

^  ^  ‘-D 
.E  c  CO 

>  ©  CD 

>  _l  _1 

t  CD  O) 

o  c  c 
OT  $  § 


11-55 


American  Institute  of  Aeronautics  and  Astronautics 

UNCLASSIFIED 


Fig.  4  -  Wind  Tunnel  Geometry  for  Computational  Fluid  Dynamics  Simulation 
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Fig.  5a  -  CFD  Simulation  of  Wind  Tunnel  Geometry  With  Windward  Jet 


Missile  Without  Thruster  Jet 
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Total  Pressure  Contours  at  Several  Axial  Locations  Showing  Jet  Effects 

Fig.  5b  -  CFD  Simulation  of  a  Missile  With  Wings,  Tail  Panels  and  Lateral  Thruster  Jet 


Missile  Without  Thruster  Jet 
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Fig.  5c  -  CFD  Simulation  of  a  Missile  With  Wings,  Tail  Panels  and  Lateral  Thruster  Jet 
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Total  Pressure  Contours  at  Various  Axial  Locations  Showing  Vortex  Formation 


Missile  With  Lateral  Windward  Thruster 
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Mach  Number  Contours  at  Several  Axial  Locations  Showing  Jet  effects 
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Fig.  7a  -  CFD  Simulation  of  Wind  Tunnel  Geometry  With  Windward  Jet 
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Fig.  7b  -  CFD  Simulation  of  Wind  Tunnel  Geometry  With  Windward  Jet 


Flow  Mach  Number-  3.94 
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Missile  Surface  Pressure  Distribution  Showing  Jet  Interaction  Effects 
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Fig.  9  -  CFD  Simulation  of  Wind  Tunnel  Geometry  With  Windward  Jet 
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Fig.  10  -  CFD  Simulation  of  Wind  Tunnel  Geometry  Without  Jet 
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Fig.  12a  -  CFD  Simulation  of  Wind  Tunnel  Geometry  With  Leeward  Jet 


Leeward  Thruster  Jet 
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Total  Pressure  Contours  at  Various  Axial  Locations  Showing  Vortex  Formation 
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Pressure  Distribution  on  Missile  Surfaces  With  Leeward  Jet  Showing  Jet  Interaction  Effects 

Fig.  13  -  CFD  Simulation  of  Wind  Tunnel  Geometry  With  Leeward  Jet 
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Fig.  14  -  CFD  Simulation  ofWind  Tunnel  Geometry  With  Leeward  Jet 


R=p/(gama*pinf) 


Fig.  15  -  CFD  Simulation  of  Wind  Tunnel  Geometry  With  Leeward  Jet 


Windward  Jet 
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Fig.  17  -  CFD  Simulation  of  Wind  Tunnel  Geometry  With  Forward  Wings 
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Fig.  19  -  CFD  Simulation  of  Wind  Tunnel  Geometry  With  Short  Wings 


Flow  Conditions:  M=3.94,  Alpha  =  20  Deg 
Windward  Jet,  Jet  Thrust=  175  Lbs 


Surface  Pressure  Differential  for  Jet  on  and  Jet  off  showing  Jet  Interaction  Effects 

Fig.  20  -  Cro  Simulation  of  Wind  Tunnel  Geometry  With  Short  Wings 
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Fig.  21b  -  CFD  Simulation  of  a  Missile  Without  Lateral  Jet 
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Fig.  21c  -  CFD  Simulation  of  a  Missile  Without  Lateral  Jet 


Missile  Rolled  45  Degree 
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Mach  Number  Distribution  at  Several  Axial  Locations  for  a  Missile 

Fig.  21d  -  CFD  Simulation  of  a  Missile  Without  Lateral  Jet 


Geometry  With  Wings,  Fins  and  Lateral  Thrusters 
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CFD  Normalized  Normal  Force  Coefficient  ( CN/CNO ) 


Missile  Geometry  with  Wings,  Fins  and  Lateral  Thruster 
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